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ABSTRACT 

We present a study of the distribution of X-ray AGN in a representative sample of 26 massive clusters at 
0.15<z<0.30, combining Chandra observations sensitive to X-ray point sources of luminosity L X ~10 42 erg s -1 
at the cluster redshift with extensive and highly complete spectroscopy of cluster members down to ~M^+2. In 
total we identify 48 X-ray AGN among the cluster members, with luminosities 2x 10 41 - 1 x 10 44 erg s -1 . Based 
on these identifications, we estimate that 0.73±0.14% of cluster galaxies brighter than Mk =-23 . 1 (M£+1.5) 
host an X-ray AGN with L^>10 42 erg s -1 . In the stacked caustic diagram that shows (vi os -<v>) / a v versus 
fproj/f5oo, the X-ray AGN appear to preferentially lie along the caustics, suggestive of an infalling population. 
They also appear to avoid the region with lowest cluster-centric radii and relative velocities (r pro j<0 Ar soo; 
|v-<v>|/cr v <0.8), which is dominated by the virialized population of galaxies accreted earliest into the clus- 
ters. The line-of-sight velocity histogram of the X-ray AGN shows a relatively flat distribution, and is incon- 
sistent with the Gaussian distribution expected for a virialized population at 98.9% confidence. Moreover the 
velocity dispersion of the 48 X-ray AGN is 1 .5 1 x that of the overall cluster population, which is consistent with 
the y/2 ratio expected by simple energetic arguments when comparing infalling versus virialized populations. 
This kinematic segregation is significant at the 4.66-cr level. When splitting the X-ray AGN sample into two 
according to X-ray or infrared (IR) luminosity, both X-ray bright (Lx>10 42 ) and IR-bright (L TIR >2x 1O 1O L ) 
sub- samples show higher velocity dispersions than their X-ray dim and IR-dim counterparts at >2a signifi- 
cance. This is consistent with the nuclear activity responsible for the X-ray and IR emission being slowly shut 
down as the host galaxies are accreted into the cluster. Overall our results provide the strongest observational 
evidence to date that X-ray AGN found in massive clusters are an infalling population, and that the cluster 
environment very effectively suppresses radiatively-efficient nuclear activity in its member galaxies. These 
results are consistent with the view that for galaxies to host an X-ray AGN they should be the central galaxy 
within their dark matter halo and have a ready supply of cold gas. 

Subject headings: galaxies: active — galaxies: clusters: general — galaxies: evolution — galaxies: stellar 
content 



1. INTRODUCTION 

In recent years observations have shown that at the heart of 
most, if not all, massive galaxies is a supermassive black hole 
(SMBH) (for a review see Ferrarese & Ford 2005). The ob- 
served tight correlations between the SMBH and galaxy bulge 
masses (e.g. Haring & Rix 2004) or stellar velocity dispersion 
(Gebhardt et al. 2000; Ferrarese & Ford 2005), imply that the 
evolution of the galaxy and the central black hole are fun- 
damentally intertwined. Silk & Rees (1998) suggested that 
these tight correlations arise naturally from the self-regulated 
growth of SMBHs through accretion of material as active 
galactic nuclei (AGN), triggered by the merger of gas-rich 
galaxies. Tidal torques produced by the merger channel large 
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amounts of gas onto the central nucleus, fuelling a powerful 
starburst and rapid black hole growth, until feedback from ra- 
diatively efficient accretion is able to drive quasar winds and 
expel the remaining gas from the galaxy (e.g. Barnes & Hern- 
quist 1992; Silk & Rees 1998; Springel et al. 2005a; Hopkins 
et al. 2008). 

There have been significant advances of late in the theoret- 
ical framework for understanding galaxy evolution by incor- 
porating AGN feedback, in particular to reproduce the global 
properties of galaxies observed in the current large scale sur- 
veys such as the SDSS (e.g. Bower et al. 2006; Croton et al. 
2006; Hopkins et al. 2006b; Schaye et al. 2010). The enor- 
mous amounts of energy produced by AGN have also been 
invoked as efficient mechanisms to quench star-formation in 
galaxies, either via extreme feedback from quasars (Springel 
et al. 2005a; Hopkins et al. 2006a,b, 2008), or via low-level, 
quasi-continuous feedback from radiatively inefficient accre- 
tion of hot gas from the X-ray emitting halos of massive 
galaxies (Best et al. 2005; Croton et al. 2006; Hardcastle et 
al. 2007). 

Given the fundamental link between SMBHs and bulge 
growth as well as the importance of AGN feedback for the 
evolution of the host galaxy, a key unresolved issue in astro- 
physics is determining the source of gas that fuels the growth 
of SMBHs and resultant nuclear activity, and what triggers 
the gas inflow. As well as the previously described merger 
paradigm, secular processes including bar-driven gas inflows 
(Kormendy & Kennicutt 2004; Jogee 2006), disk-instabilities, 
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stochastic collisions with molecular clouds (Hopkins & Hern- 
quist 2006), and stellar winds from evolved stars (Ciotti & 
Ostriker 2007) have all been proposed as means to supply gas 
onto SMBHs and trigger their activity. 

A primary complication for understanding the nature of 
AGN and their impact on galaxy evolution is the observed 
zoo of AGN populations and classes. These are variously 
identified via optical spectroscopy or characteristic emission 
at X-ray, radio or infrared wavelengths, and each one shows 
a diverse range of Eddington ratios or host galaxy prop- 
erties. In principle, X-ray emission provides the cleanest 
and most reliable approach to selecting AGN and in partic- 
ular identifying those galaxies undergoing significant black 
hole growth. The hard X-ray emission is directly associ- 
ated with accretion close in to the black hole (10-100 grav- 
itational radii). It is produced in the hot corona that sur- 
rounds the black hole by the inverse Compton scattering of 
ultraviolet photons emitted by the accretion disk. Except for 
strongly absorbed AGN (N#>10 23 cm -2 ) which are rare at 
low-redshifts (Tozzi et al. 2006; Burlon et al. 2011), X-ray 
selection is nearly independent of obscuration. In contrast the 
selection of AGN at UV/optical wavelengths is strongly af- 
fected by absorption and orientation effects, as well as being 
prone to contamination by stellar light. Moreover, while X- 
ray and infrared- selected AGN typically have Eddington ra- 
tios of ~1-10% (Hickox et al. 2009), associating them with 
rapid black hole growth, radio-loud AGN and low-luminosity 
optically- selected AGN (LINERs) are characterized by much 
lower Eddington ratios (< 10~ 3 ). 

The environments, galaxy hosts and clustering of the di- 
verse classes of AGN can provide important clues to under- 
standing both the accretion processes powering the AGN (e.g. 
whether it is fuelled by the accretion of hot or cold gas), and 
its subsequent impact on the galaxy host in terms of building 
up its bulge or quenching its star-formation. This has moti- 
vated numerous authors who have examined the AGN pop- 
ulations in massive galaxy clusters (e.g. Dressier, Thompson 
& Schectman 1985; Martini et al. 2006, 2007, 2009) where, 
for the virialized galaxy population at least, only hot gas is 
available for accretion onto the SMBHs. 

The distribution of AGN in galaxy clusters provides a fun- 
damental test for those accretion processes which require a 
ready supply of cold gas in the host galaxy, such as secu- 
lar bar/disk-instabilities (Hopkins & Hernquist 2006) or the 
merger of two gas-rich galaxies. Given the much lower frac- 
tion of galaxies with substantial reservoirs of Hi gas found in 
cluster cores compared to the field (e.g. Giovanelli & Haynes 
1985; Cortese et al. 2011), nuclear activity among the viri- 
alized cluster population should also be reduced. However, 
the molecular gas contents of late-type cluster galaxies ap- 
pear normal (Boselli & Gavazzi 2006), while Young et al. 
(2011) and Haines et al. (2011) show that at least some clus- 
ter SOs are able to retain their molecular gas and dust con- 
tents through several Gyr while virialized within the cluster. 
As disk-instabilities are described by the Toomre criterion, it 
may be sufficient to simply reduce the gas surface density of 
cluster galaxies to prevent secular triggering of AGN, by mak- 
ing the remaining gas stable against collapse. In the case of 
the merger paradigm, nuclear activity should be even more 
strongly suppressed within rich clusters. Here the encounter 
velocities of galaxies are much greater than their internal ve- 
locity dispersions, preventing their coalescence, in spite of the 
high galaxy densities (Aarseth & Fall 1980). Gas-rich merg- 
ers should instead be most frequent in galaxy groups (Hop- 



kins et al. 2008) and the cluster outskirts, where many galax- 
ies (including gas-rich ones), under the influence of the clus- 
ter's tidal field, are part of a convergent flow resulting in en- 
hanced interactions between neighbors (van de Weygaert & 
Babul 1994). Since these galaxies are falling into the clus- 
ter for the first time, they have yet to have their gas contents 
be affected by their passage through the dense intra-cluster 
medium (ICM). Galaxy harassment due to frequent high- 
speed fly-by interactions has also been proposed as a means 
of triggering nuclear activity, by driving dynamical instabili- 
ties that efficiently channel gas onto the SMBHs (Moore et al. 
1996), although again this requires the host galaxy to contain 
a gas reservoir. 

The trends predicted by associating X-ray AGN to gas-rich 
galaxies could be diluted however, or even reversed, by the 
tendency of X-ray AGN (above a given L x ) to be hosted by the 
most luminous galaxies (Sivakoff et al. 2008; Tasse, Rottering 
& Best 2011), a population which is most centrally concen- 
trated within clusters (Lin, Mohr & Stanford 2004; Thomas 
& Katgert 2006). Both Sivakoff et al. (2008) and Haggard et 
al. (2010) show the X-ray AGN fraction to increase by an or- 
der of magnitude from low-mass populations (7W^1O 1O M ) 
to the most massive galaxies (M>10 U M Q ) in both cluster 
and field populations. These trends act as physical selection 
effect, as the most luminous galaxies are more likely bulge- 
dominated and hence have higher black hole masses. Thus 
for the same accretion rate relative to the Eddington limit, the 
higher mass galaxy, will likely have a higher X-ray luminos- 
ity, and be more likely detected above a fixed L x limit. 

If nuclear activity is not dependent on a cold gas supply, but 
is being instead fuelled by the accretion of hot gas from the 
halos of massive galaxies (Croton et al. 2006) or recycled gas 
from evolved stars (Ciotti & Ostriker 2007) then we would ex- 
pect nuclear activity to be most prevalent among the massive 
elliptical galaxies that dominate cluster cores. In this case the 
AGN fraction should be higher in clusters than the field. 

The earliest works identified AGN via their optical emis- 
sion lines (Gisler 1978; Dressier, Thompson & Schectman 
1985), finding them to be much rarer among cluster mem- 
bers than in field samples. More recent studies based on the 
SDSS have obtained rather conflicting results, with some con- 
firming the previous decline in the AGN fraction with galaxy 
density (Kauffmann et al. 2004; Mahajan et al. 2010; Hwang 
et al. 2012), while others find no apparent trend with environ- 
ment (Miller et al. 2003; Haines et al. 2007). von der Linden 
et al. (2010) find that while overall the fraction of galaxies 
which host a powerful optical AGN (with L[Olll]>10 7 L o ) 
decreases towards the cores of clusters, when considering 
just star-forming galaxies, the AGN fraction remains constant 
with cluster-centric radius. They suggest that the decline in 
the AGN fraction is almost exclusively an effect of the in- 
creasing number of red sequence galaxies, a population which 
hosts very few strong AGN. 

The launches of the Chandra and XMM X-ray Observa- 
tories have opened up an efficient way of identifying X-ray 
AGN in and around massive clusters from the same images 
used to probe the fundamental properties of the clusters them- 
selves (e.g. masses, 7)cm). Statistical analyses of X-ray point 
sources in cluster fields have reported overdensities of X-ray 
sources with respect to non-cluster fields (e.g. Cappi et al. 
2001 ; Molnar et al. 2002; Ruderman & Ebeling 2005; Branch- 
esi et al. 2007). Gilmour et al. (2009) found an excess of ~ 1 .5 
X-ray point sources per cluster within 1 Mpc from a sample 
of 148 clusters at 0.1<z<0.9. They found the radial distri- 
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bution of these excess sources consistent with a flat radial 
distribution within 1 Mpc, although they could also be con- 
sistent with being drawn from the general cluster galaxy pop- 
ulation. Ruderman & Ebeling (2005) found instead a central 
spike within 0.5 Mpc followed by a secondary broad peak at 
the virial radius for relaxed clusters, with flatter spatial distri- 
butions for disturbed systems. Koulouridis & Plionis (2010) 
instead found that the overdensity of X-ray point sources in 
16 clusters at 0.07<z<0.28 to be a factor ~4 less than that 
of bright optical galaxies, and concluded that the triggering of 
luminous (Lx>10 42 ergs -1 ) X-ray AGN to be strongly sup- 
pressed in rich clusters. 

Martini et al. (2006) performed a redshift survey of X- 
ray point sources in 8 clusters at 0.06<z<0.31, identi- 
fying 30 X-ray AGN with L x >10 41 ergs -1 and M*<-20. 
Six of these clusters were X-ray luminous systems with 
L x >4xl0 44 ergs -1 , including Abell 1689. The other two 
were the pair of low-mass merging clusters Abell 3125 and 
Abell 3128 at z=0.06. Abell 3125 in fact does not have a dif- 
fuse ICM, while the ICM for Abell 3128 is double peaked. 
The resulting X-ray AGN fraction of ~5% among cluster 
galaxies brighter than M R =-20, was much higher than pre- 
viously expected and remarkably similar to those found in the 
field (Martini et al. 2007; Lehmer et al. 2007). This suggested 
little or no environmental dependence and that cluster galaxies 
could retain significant reservoirs of cold gas near their cen- 
tral SMBHs. The X-ray AGN fraction was highest (-12%) 
in the two low-mass merging systems Abell 3125/3128. Only 
a small fraction of these X-ray AGN would have been clas- 
sified as such from their optical spectra, and so the differ- 
ence in the observed environmental trends seen in optically- 
selected QSOs in which they avoid cluster cores (e.g. Dressier, 
Thompson & Schectman 1985; Sochting et al. 2004), was put 
down to biases against detection of low-luminosity AGN due 
to host galaxy dilution and obscuration. In the dynamical 
analysis of the same sample, both the stacked velocity and ra- 
dial distributions of the X-ray AGN were found to be entirely 
consistent with being drawn from the overall cluster galaxy 
population (Martini et al. 2007). The most luminous X-ray 
AGN with L X >10 42 ergs -1 were however found to be more 
centrally concentrated than inactive cluster galaxies. This 
study had little coverage of the cluster outskirts, with 90% 
of the X-ray AGN within 0.5 r2oo- However, Sivakoff et al. 
(2008) using mosaics of Chandra images covering the virial- 
ized regions of Abell 85 and Abell 754 (two merging clusters 
at z=0.055 with Lx~5 x 10 44 erg s -1 ), obtained similar trends, 
with again no differences in the spatial distributions (both ra- 
dially or in velocity) of X-ray AGN (with L^>10 41 ergs -1 ) 
and non-active cluster galaxies. Gilmour et al. (2007) instead 
found that the 1 1 X-ray AGN within the supercluster A901/2 
at z~0.17 lie predominately in intermediate density regions, 
avoiding the cluster cores, prefering instead the cluster out- 
skirts or group-like environments containing more blue galax- 
ies than on average. 

From these analyses it remains unclear whether the X- 
ray AGN observed in galaxy clusters should be considered 
to be virialized or an infalling population. Given the small 
numbers of clusters studied in detail with extensive spec- 
troscopy, it is unlikely that the previous three samples can be 
regarded as representative of the general massive cluster pop- 
ulation. The samples of Martini et al. (2006), Sivakoff et al. 
(2008) and Gilmour et al. (2007) all appear overrepresented 
by merging systems. What is lacking is a statistical sample 
of spectroscopically-confirmed X-ray AGN from a large, rep- 



resentative set of massive clusters at a specific epoch. The 
aim of this present work is to produce such a sample, tak- 
ing advantage of the unique multi- wavelength dataset recently 
assembled within LoCuSS (Local Cluster Substructure Sur- 
vey) 8 . 

LoCuSS is a multi- wavelength survey of X-ray luminous 
galaxy clusters at 0.15<z<0.3 drawn from the ROSAT All 
Sky Survey cluster catalogs (Ebeling et al. 1998, 2000; 
Bohringer et al. 2004). LoCuSS comprises three broad 
inter-related themes: (i) calibration of mass-observable scal- 
ing relations for cluster cosmology, (ii) physics of clus- 
ter cores, and (iii) the physics of galaxy transformation in 
group/cluster environments. This is the latest in a series of 
papers in the third "galaxy evolution" theme. The science 
goals and design of this theme are described in detail by Smith 
et al. (2010b). In summary, we have obtained panoramic 
multi- wavelength data (FUV-FIR) on a morphologically un- 
biased sample of 30 clusters whose X-ray luminosities span 
3.79xl0 44 <L x (0.1-2.4keV)<2.28xl0 45 ergs -1 . These data 
span fields of >25 / x25 / - i.e. out to ~1.5-2 cluster virial 
radii (Haines et al. 2009a; Pereira et al. 2010; Smith et al. 
2010b). These 30 clusters were selected from the parent sam- 
ple simply on the basis of being observable by Subaru on the 
nights allocated to us (Okabe et al. 2010), and should there- 
fore not suffer any gross biases towards (for example) cool 
core cluster, merging clusters etc. Indeed Okabe et al. (2010) 
show that the sample is statistically indistinguishable from a 
volume-limited sample. 

Twenty-six of the 30 clusters have been observed with 
Chandra to depths suitable for identifying X-ray AGN within 
the clusters. We have combined these X-ray data with our ex- 
tensive spectroscopy (^200-400 members per cluster) from 
MMT/Hectospec, to derive a highly complete census of the 
X-ray AGN population within this statistical sample of local 
clusters. We have performed a dynamical analysis of these 
cluster X-ray AGN in order to understand whether they are 
hosted primarily by galaxies still infalling into the cluster (and 
which may have yet to encounter the dense ICM) or rather 
galaxies which were accreted at early epochs and are now viri- 
alized. This analysis is aided by comparison to predictions 
from the Millennium cosmological simulation of the distri- 
butions of galaxies in the stacked caustic diagrams coded ac- 
cording to the epoch at which they were accreted into the clus- 
ter. We also present the X-ray luminosity function of cluster 
AGN and examine the global optical and infrared properties 
of the galaxies hosting X-ray AGN in clusters. 

In §2 we summarize the data used in this paper, while our 
use of cosmological simulations is described in §3. The main 
results are then presented in §4, followed by a discussion and 
summary in §5 and §6. Throughout we assume ^m=0.3, 
f^A=0.7 and H =70kms -1 Mpc -1 , for consistency with our 
previous work (Haines et al. 2009a). 

2. DATA 

This sample of 26 clusters was chosen because the clusters 
all have publicly available Chandra ACIS data with exposure 
times of ~10 ksec or greater, sufficient to detect X-ray sources 
as faint as L^^IO 42 erg s -1 at the cluster redshift, and to de- 
rive robust measures of rsoo from the extended X-ray emission 
from the ICM. Of the original sample of 30 clusters, A291 and 
A2345 were excluded due to lack of Chandra imaging, while 
it was not possible to determine rsoo values for A689 or Z348. 

8 www.sr.bham.ac.uk/locuss/ 
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For each of the clusters we have panoramic near-IR imaging 
to K~19, J~21, obtained with either WFCAM on the 3.8-m 
United Kingdom Infrared Telescope (UKIRT) which covers 
52' x 52' regions centered on each cluster (23 out of 26 clus- 
ters), or NEWFIRM on the 4.0-m telescope at Kitt Peak which 
provides a field of view of 27 / x27 / for the remaining 3 sys- 
tems (Haines et al. 2009a). For each cluster we also have 
deep Subaru optical (V,/) imaging covering 34 / x27 / , used 
by Okabe et al. (2010) to perform weak lensing mass esti- 
mates. Additionally, ugriz photometry to r~22 is available 
for 22 of the 26 clusters, taken from the eighth data release 
from the Sloan Digital Sky Survey (SDSS-DR8; Aihara et al. 
2011). Each cluster was observed across a 25 / x25 / field of 
view at 24/im with MIPS (Rieke et al. 2004) on board the 
Spitzer Space Telescope 9 (Werner et al. 2004) (PID: 40872; 
PI: G. P. Smith), reaching typical 90% completeness limits 
of 400/iJy, sensitive to obscured star formation within clus- 
ter galaxies to levels of 3 M yr -1 in our most distant sys- 
tems. Herschel PACS/SPIRE imaging from our LoCuSS Her- 
schel Key Programme provides 100-500/im photometry for 
matched 25'x25' fields (Smith et al. 2010b), reaching com- 
parable sensitivities in terms of obscured SFRs. 

2.1. MMT/Hecto spec spectroscopy 

We have recently completed ACReS (the Arizona Cluster 
Redshift Survey; Pereira et al. in preparation) a long-term 
programme to observe our sample of 30 galaxy clusters with 
MMT/Hectospec. Hectospec is a 300-fiber multi-object spec- 
trograph with a circular field of view of 1 ° diameter (Fabri- 
cant et al. 2005) on the 6.5m MMT telescope at Mount Hop- 
kins, Arizona. The details of the observations, targeting strat- 
egy and data reduction are presented elsewhere (Haines et al. 
2009a,b; 2010; Pereira et al. in preparation). 

Briefly, probable cluster galaxies are targeted according to 
their AT-band absolute magnitude and J-K color, prioritizing 
those galaxies detected with Spitzer/MIPS at 24/im. The tar- 
geting is based on the empirical observation that galaxies of 
a particular redshift lie along a single narrow J-K/K color- 
magnitude (C-M) relation, which evolves redward monoton- 
ically with redshift to z~0.5 (Haines et al. 2009b). The NIR 
colors of galaxies are relatively insensitive to star-formation 
history and dust extinction. This effectively allows the cre- 
ation of a stellar mass-limited sample within a narrow red- 
shift slice centered on the cluster, with no preference towards 
passive ("red sequence") or star-forming ("blue cloud") sub- 
populations, simply by selecting galaxies within a color slice 
of width 0.3-0.4 mag enclosing the observed J-K/K C-M re- 
lation for that cluster (Haines et al. 2009a). For each clus- 
ter, galaxies brighter than M£ +2.0 were targeted (falling to 
+1.5 for particularly rich systems at higher redshifts) over 
the full UKIRT/WFCAM field, reaching overall completeness 
levels (as measured within r 2 oo) of 70% for the M^-selected 
sample, increasing to 96.4% (5245/5441) for those galaxies 
detected with Spitzer. 

The selection of galaxy targets was made without any prior 
knowledge of whether they were detected with Chandra or 
not, and so should not have any bias with respect to X-ray 
AGN. However, X-ray AGN are more likely to be MIR-bright 
than normal galaxies of the same stellar mass (Treister et al. 
2006), and so are more likely to have been targeted. This 

9 This work is based in part on observations made with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology under a contract with NASA (contract 1407) 



Table 1 

Summary of Chandra observations. 
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2 
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5 
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1 
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RXJ2129.6+0005 


0.2337 


1.227 
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1 
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was accounted for explicitly by weighting each galaxy ac- 
cording to the probability that it was targeted for spectroscopy 
based on its Mk, 24/im flux and spatial location within the 
cluster. Moreover, our spectroscopic survey should have no 
morphological bias against unresolved extra-galactic objects 
such as compact galaxies or quasars. This is because we tar- 
geted sources based solely upon their near-IR colors or 24/im 
emission, irrespective of whether they are resolved or not in 
our near-IR images. Stars show much bluer near-IR colors 
(J-K<\) than galaxies or quasars at the redshift of interest 
(see Figs 2, 3 of Haines et al. 2009b), and were identified and 
excluded as targets for spectroscopy only if they have J-K< 1 
and are unresolved in our AT-band data. 

To date, ACReS has required the equivalent of 13 full nights 
of observations since December 2008, producing ^30000 
spectra, of which ~10 000 have been identified as being clus- 
ter members, the largest sample to date from a redshift survey 
targeting galaxy clusters. We have obtained 3-6 fiber config- 
urations per cluster resulting in typically redshifts for 1 00- 
500 cluster members, the number depending primarily on the 
richness and/or compactness of the cluster. The 270 line grat- 
ing was used, providing coverage over the wavelength range 
3 65 0-9200 A at 6. 2 A resolution. For each cluster, we plot- 
ted redshift against cluster-centric radius, identifying cluster 
members as lying within the general "trumpet" -shaped caus- 
tic profile expected for galaxies infalling and orbiting within 
a massive structure (see §3). 

2.2. Chandra X- ray imag ing 
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Fig. 1. — Broadband (0.3-7 keV) X-ray counts of detected point sources 
with known redshifts as a function of projected clustercentric radius, reveal- 
ing the impact of the extended X-ray emission from the ICM on point source 
detection. Solid (open) symbols indicate X-ray point sources associated to 
cluster members (non-members). The horizontal dashed line indicates the 
minimum requirement for detection of 6 counts. 
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Fig. 2. — Positional offsets (A RA, A Dec) between the X-ray point source 
and the ^-band counterpart for all counterparts with known redshifts. Solid 
symbols indicate cluster members. The size of the symbol scales with the 
number of X-ray counts. 



The 26 clusters have available Chandra data, both from 
the archive and drawn from our own Chandra Cycle 10 pro- 
gramme (PID: 10800565). The observations of 20/26 clusters 
were made with the I mode of the Advanced Camera for Imag- 
ing Spectroscopy (ACIS-I), the remaining six being observed 
with ACIS-S. The exposure times were typically 20ksec, but 
range between 9-120ksec. The larger ACIS-I field of view 
is 16.9' x 16.9', slightly smaller than our SpitzerMWS fields, 
meaning that all sources detected by Chandra are also covered 
by our 24/nm imaging. The ACIS-I instrument is made up of 
4 CCDs in a 2x2 grid. The cluster core is always centered 
within one of these CCDs to minimize the impact of the gaps 
between the CCDs on measurements of the ICM, and hence 
the spatial coverage of the clusters is rather asymmetric. The 
ACIS-S imaging provides a smaller field of view of 8. 3 7 x 8.3' 
from a single CCD, centered on the cluster core. 

The deprojected dark matter densities, gas densities and gas 
temperature profiles were derived by fitting the phenomeno- 
logical cluster models of Ascasibar & Diego (2008) to a se- 
ries of annular spectra extracted for each cluster (Sanderson 
& Ponman 2010). The best-fitting cluster models were then 
used to estimate rsoo, the radius enclosing a mean overdensity 
of 500 with respect to the critical density of the Universe at 
the cluster redshift (Sanderson et al. 2009). A summary of the 
clusters observed, redshifts and rsoo values, Chandra observa- 
tion ID and ACIS instrument used and total exposure times 
after cleaning is shown in Table 1 . Details of the reductions 
and analysis can be found in Sanderson et al. (2009). 

To detect X-ray point- sources that are potential X-ray AGN, 
the wavelet-detection algorithm CIAO WAVDETECT was ap- 
plied, requiring a minimum of six counts in the broad en- 
ergy (0.3-7 keV) range to be considered in this analysis. The 
observed fluxes in this band were derived assuming a T=1.7 
power-law spectrum with Galactic absorption, following Ren- 
ter et al. (2005). We then calculated the rest- frame luminosity 
for all sources with redshifts, assuming k-corrections of the 
form(l+z) r ~ 2 . 

For our typical exposure time of 20 ksec, our survey sensi- 



tivity limit of six counts corresponds to an X-ray luminosity 
of Lx=3.5xl0 41 ergs _1 at z=0.20. Considering the range of 
combinations of exposure times and redshifts for our Chan- 
dra cluster survey, we find our on-axis survey limits are at 
or below L x =l .0 x 10 42 erg s -1 at the cluster redshift for all 26 
systems (column 6 of Table 1). We expect this sensitivity to 
reduce by no more than 10-20% for sources with the largest 
off-axis angles (~10') 10 

The presence of extended emission from the ICM affects 
the ability of WAVDETECT to detect faint X-ray point sources 
near the cluster core due to increased photon noise. This effect 
is demonstrated in Fig. 1 which shows the number of broad- 
band X-ray counts as a function of cluster-centric radius for 
each point source matched to a galaxy counterpart with known 
redshift. Within 0.4 r 50 o (0.2 r 50 o) WAVDETECT doesn't detect 
any point sources with fewer than 10 (20) counts, indicating a 
~2-3 x reduction in sensitivity to X-ray point sources in the 
cluster cores. The impact of the ICM background on our abil- 
ity to detect X-ray AGN varies significantly from cluster to 
cluster, being significantly worse for the cool-core clusters for 
which the X-ray emission is strongly concentrated in the cores 
(and which tend to have deeper images) than their non cool- 
core counterparts. We confirm that for 20/26 clusters we re- 
main complete to at least L x =l .0 x 10 42 erg s -1 even in the core 
regions with the strongest ICM emission. We note also that 
all five cluster X-ray AGN identified within 0.2r5oo have X- 
ray luminosities of 0.87-2.2 x 10 42 erg s -1 . For the remaining 
six clusters A611, A1835, A2390, RXJ 2129.6+0005, Z2089 
and Z7160 their concentrated X-ray emission related to their 
strong cool cores does means that we are locally (only within 
<0.1 r5oo) insensitive to X-ray AGN at this luminosity level. 

We identified possible NIR counterparts to the X-ray point 
sources, matching objects within a generous limit of 4 arcsec. 
Given the positional accuracy of the X-ray centroids in Chan- 
dra images, true counterparts are expected to be within 
2 arcsec, except in the cases of large, extended galaxies for 
which the optical centroid may not coincide well with the ac- 

10 Chandra Proposers' Guide: http://cxc.harvard.edu/proposer/. 
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tive nucleus. Each possible match was then visually checked 
in our deep Subaru optical images, and in the case of multi- 
ple possible counterparts, the nearest match taken, unless the 
other is infrared-bright and within 2 arcsec. Figure 2 shows 
the positional offsets (A RA, A Dec) for all the optical coun- 
terparts to X-ray point sources with known redshifts. The vast 
majority have counterparts within 1.5 arcsec, including 46 of 
48 cluster members (shown as solid symbols), confirming that 
the associations are reliable. 

Over the full sample of 26 clusters we have redshifts for 
151 out of the 183 galaxies detected with Chandra that would 
be brighter than M£+2 at the cluster redshift, and hence could 
have been targeted for spectroscopy within ACReS. This gives 
us a spectroscopic completeness of 83% for the X-ray AGN 
subsample. We note however that only three of the remain- 
ing 32 X-ray sources would be identified as targets for spec- 
troscopy within the ACReS selection criteria (§ 2.1). The 
rest are classified as being too red in J-K to be cluster mem- 
bers (being on average ~0.2mag redder than our J-K color 
cut and ~0.4mag redder than the cluster C-M relation), and 
instead are likely background galaxies. This is supported 
by them being on average ~ 1.3 mag fainter in the AT-band 
than our spectroscopic cluster X-ray AGN sample (which has 
a median magnitude of ^=15.2). None of these 32 X-ray 
sources lacking redshifts would be brighter than L* at the clus- 
ter redshift, whereas this is the median ^-band luminosity of 
our confirmed cluster X-ray AGN sample. Thus, our effective 
spectroscopic completeness is >95% for X-ray point sources 
hosted by galaxies brighter than M£ +2 whose near-IR colors 
are consistent with being at the cluster redshift. In terms of 
spatial coverage, we find that our Chandra imaging covers 
90% of cluster members within rsoo or 75% of those within 
?*200 over the full sample. 

We have spectroscopically identified a total of 48 member 
galaxies with X-ray emission in 26 clusters. For those host 
galaxies detected with Spitzer at 24/im, we estimated their star 
formation rates (SFRs) and total infrared luminosities (L TIR ) 
by fitting the infrared-SED models of Rieke et al. (2009) to the 
observed 24/im fluxes. Table 2 lists their positions, redshifts 
and X-ray luminosities along with a summary of the key prop- 
erties (M K , SFR, morphological and spectral classes) of their 
host galaxies. BCGs (defined here as the Brightest Cluster 
Galaxy in the cluster core) have been explicitly excluded due 
to their unique star formation histories and evolutions (Lin & 
Mohr 2004), and the direct link between BCG activity and the 
presence of cooling flows within clusters (Edge 1991; Bildfell 
et al. 2008; Smith et al. 2010a). It is also unclear whether the 
X-ray emission comes from the galaxy itself or the ICM, par- 
ticularly in the case of cooling flow clusters. 

3. THE MILLENNIUM SIMULATION AND THE 
DISTRIBUTION OF GALAXIES IN THE CAUSTIC 
DIAGRAM AS A FUNCTION OF ACCRETION EPOCH 

The location of infalling and virialized populations in and 
around massive clusters are well separated in radial phase- 
space diagrams (r, v ra diai\ Mamon et al. 2004; Dunner et al. 
2007; Mahajan et al. 2011). This is not the case for the ob- 
servable counterpart of projected cluster-centric radius (r pro j) 
versus line-of-sight (LOS) velocity (vi os -<v>) relative to the 
cluster redshift. However, via the use of cosmological simu- 
lations containing dozens of massive clusters similar to those 
in LoCuSS, populated by galaxies based on semi-analytic 
models, it is possible to produce simulated caustic diagrams 
for which the accretion histories of the member galaxies are 



known. For this purpose, 20x20x40/T 3 Mpc 3 regions cen- 
tered on the 30 most massive clusters were extracted from 
the Millennium Simulation (Springel et al. 2005a). These 
simulations cover a (500/z _1 Mpc) 3 volume, producing dark 
matter (DM) halo and galaxy catalogs based on the semi- 
analytic models (GALFORM) of Bower et al. (2006) for which 
positions, peculiar velocities, absolute magnitudes and halo 
masses are all provided at 63 snapshots to z=0. This allows 
the orbit of each galaxy with respect to the cluster to be fol- 
lowed from formation up until the present day, from which the 
epoch at which it is accreted (here defined as when it passed 
within r5oo(z) for the first time) can be determined. The virial 
masses of these 30 cluster halos at z=0.21 cover the range 
2.8-21.6x 1O 14 /T 1 M , consistent with the mass range of 2.9- 
14.7xlO 14 /z _1 M for our cluster sample based on the weak 
lensing mass estimates of Okabe et al. (2010). 

Knowing the relative position and peculiar velocity of each 
galaxy with respect to its host cluster it is easy to reproduce 
observations of that cluster in the form of projected cluster- 
centric radii and LOS velocity. The 30 systems are then 
stacked, scaling each by the cluster radius (rsoo) and veloc- 
ity dispersion, a v . This stacking over many clusters is not re- 
quired simply to match our observed sample, or even to con- 
tain sufficient galaxies to derive robust statistical properties 
(although this is important), but most importantly to even out 
all of the large line-of-sight variations due to the presence of 
infalling structures or filaments for individual cluster lines of 
sight and times of observation. This applies both for the simu- 
lations and the observations, motivating our use of such large 
cluster samples. 

Figure 3 shows the stacked caustic plot for galaxies orbit- 
ing the 30 most massive clusters in the Millennium simula- 
tion, as would be observed at z=0.21. The x-axis shows the 
cluster-centric radius projected along the z-direction of the 
simulation, scaled by r5 o(z), while the y-axis shows the LOS 
velocity (along the z-axis of the simulation) of each galaxy 
scaled by the velocity dispersion of galaxies within that clus- 
ter. Each symbol indicates a M^<-23.1 galaxy colored ac- 
cording to the epoch when it was accreted by the cluster. The 
overall distribution of these galaxies forms the characteris- 
tic "trumpef'-shaped caustic profile expected for galaxies in- 
falling and subsequently orbiting within a massive virialized 
structure (Regos & Geller 1989; Dunner et al. 2007). 

The distribution of galaxies in Figure 3 suggests that 
broadly speaking, the galaxies in clusters can be thought of 
as belonging to two broad categories. The first population are 
those galaxies accreted at an early epoch, and are shown in 
Figure 3 as red symbols. These are spatially localized in the 
cluster core and their LOS velocities are typically <lcr. We 
identify this population with those galaxies that either formed 
locally or were accreted when the cluster's core was being as- 
sembled, with their low velocities reflecting the fact that the 
system they fell into was only a fraction of the present-day 
mass of the cluster and possibly further slowed by dynamical 
friction. The second population of galaxies are those accreted 
after the formation of the core. Some of these galaxies are 
falling in for the very first time. These galaxies span the gamut 
from those that have yet to cross rsoo (blue symbols), those 
that are close to pericenter, and through to "backsplash galax- 
ies" (green symbols, primarily) that have completed their first 
pericenter and are now outward bound (Mamon et al. 2004; 
Pimbblet 2011). The maximum extent of this population in 
the vertical axis of Figure 3 is a function of both the actual 
velocities of the galaxies and a geometric projection factor. 
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Fig. 3. — Stacked phase-space diagram (r pro j/rsoo vs (v/ os -<v>)/cr v ) for the 30 most massive clusters in the Millennium simulation at z=0.21, the snapshot 
closest to the mean redshift of our cluster sample. Each point indicates a M#<-23.1 galaxy from the Bower et al. (2006) semi-analytic model catalog, colored 
according to when it was accreted onto the cluster, those being accreted earliest indicated by red symbols. Here accretion epoch is defined as the snapshot at 
which the galaxy passes within r^oo for the first time. Galaxies yet to pass through r^oo by z=0.21 are indicated in mid-blue and dominate at large cluster-centric 
radii (r>2r5oo) and along the caustics |v— <v> |/cr v >l). 



At small projected radii, galaxies with the highest LOS veloc- 
ities are physically located deep inside the cluster, their high 
velocities being the result of being accelerated all the way in. 

Figure 3 demonstrates that we can statistically asso- 
ciate those galaxies lying along the caustics as mostly in- 
falling (blue symbols), and those with low LOS velocities 
and cluster-centric radii as the most likely to be virialized 
(red/orange symbols). It reveals the potential of using the 
caustic diagram to statistically identify the primary accretion 
epoch(s) of observed cluster galaxy populations based upon 
their distribution in the plot. 

4. RESULTS 

In total we have identified 48 X-ray point sources among 
member galaxies in the 26 clusters studied, corresponding to 
between and 6 X-ray sources per cluster (Table 1). Of these, 
24 have both L x >10 42 ergs- 1 andM^<-23.1 (K*+1.5). From 
our spectroscopic survey of other cluster members within 
the Chandra images, we identify in total 2702 M^<-23.1 
galaxies, giving us an estimate of fraction of massive cluster 
galaxies (M^<-23.1) hosting X-ray AGN (L x >10 42 ergs _1 ) 
of 0.7 3 ±0.1 4 percent, once we account for the fact that the 
spectroscopic completeness of the X-ray AGN subsample is 
higher than that for the remaining inactive cluster popula- 
tion. This is consistent with Martini et al. (2007) who found 
that 0.87 ±0.25 percent of M R <-20 galaxies in eight clusters 



host AGN with L^>10 42 erg s" 1 . Although their survey is R- 
band selected rather than our ^-band selection, both reach 
r^M* +1.5 and so should be directly comparable. Interest- 
ingly, our value is ~40% lower than the 1.19±0.11 percent 
obtained by Haggard et al. (2010) for M R <-20 field galaxies 
at 0.05<z<0.31 from their analysis of 323 archive Chandra 
images covered by SDSS DR5 imaging and spectroscopy. 

Thirty-two of the 48 X-ray AGN are detected with Spitzer 
and nine of these are identified as LIRGs (Luminous InfraRed 
Galaxies; L TIR >10 n L (I) , where TIR is the Total InfraRed lu- 
minosity integrated over 8-1000/im). Conversely, 6 out of 
the 10 most luminous cluster galaxies at 24 fim (covered by 
the Chandra data) are found to be X-ray AGN. We note 
that this does not mean that these IR-bright AGN would be 
identified as AGN via the infrared selection of Stern et al. 
(2005), while the 24/im emission may be due to star forma- 
tion rather than nuclear activity. We do note however that for 
eight of the nine LIRGs, a comparison to our available Her- 
schel photometry reveals that they have much flatter SEDs 
((/ioo//24)<10) than found for star-forming galaxies at these 
redshifts ((/ioo//24)~25; Smith et al. 2010b; Pereira et al. 
2010), and more consistent with the flatter power-law SEDs of 
Polletta et al. (2007) in which the infrared emission is mostly 
due to dust heated by nuclear activity. For the remaining 
24/im-detected X-ray AGN, those objects detected by Her- 
schel mostly have far-IR SEDs consistent with star-forming 
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Fig. 4. — Broadband X-ray luminosity, Lx, versus the near-infrared abso- 
lute magnitude, Mk, for X-ray detected galaxies at 0.15<z<0.30. X-ray loud 
cluster members detected (not-detected) at 24/im are indicated by magenta 
(black) symbols, while field galaxies are shown as open blue (black) sym- 
bols. Larger symbols indicate galaxies identified as LIRGs from their 24 ^m 
fluxes. The expected levels of X-ray emission from low-mass X-ray binaries 
(dot-dot-dot-dashed line; Kim & Fabbiano 2004) or thermal emission from 
the diffuse gas halo (dot-dashed line; Sun et al. 2007) are shown. 
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Fig. 5. — The broadband (0.3-7 keV) X-ray luminosity function of X-ray 
AGN belonging to the 26 clusters. The right-hand axis shows numbers of 
X-ray AGN observed per 0.5 dex bin in luminosity, while the left-hand axis 
shows the averaged number of X-ray AGN per cluster within r2oo (including 
a 25% correction for those galaxies within r2oo not covered by Chandra data) 
per dex in Lx. The dashed line shows the z=0.3 field X-ray luminosity func- 
tion of Hasinger et al. (2005) matched to have the same volume as our cluster 
sample (3.1 x 10 4 Mpc 3 ). It has been evolved in space-density to z=0.225 and 
luminosity-corrected from 0.5-2 keV to 0.3-7 keV assuming a T=1.7 power- 
law spectrum. The solid curve shows the same luminosity function multiplied 
by a factor 7.5 to match the cluster LF. 



galaxies. 

At X-ray luminosities L^<10 42 erg s -1 there are four poten- 
tial sources of X-ray emission in galaxies: low-mass X-ray 
binaries (LMXBs) which are sensitive to the total stellar mass 
of the galaxy (Kim & Fabbiano 2004); high-mass X-ray bina- 
ries (HMXBs) which are sensitive to the recent star formation 
activity (Mineo, Gilfanov & Sunyaev 2012); thermal emis- 
sion from a hot gas halo; and an AGN. It it thus important at 
this point to quantify the possible contribution to the X-ray 
emission of our Chandra point sources from the former three 
components. Mineo, Gilfanov & Sunyaev (2012) show that 
HMXBs are a good tracer of recent star formation, and their 
collective luminosity scale linearly with the SFR, L x (0.5- 
8.0keV)-2.5 x 10 39 SFR(M yr _1 ) for SFRs in the range 0.1- 
1000 M yr -1 . Based on comparison of the 24/im-based SFR 
estimates to the observed X-ray luminosities for each point 
source, HMXBs contribute a maximum of 10-20% of the X- 
ray luminosity, and in most cases less than 1%. To examine 
the possible contribution to the X-ray emission from LMXBs 
or the hot gaseous coronae, which should both scale with stel- 
lar mass, Figure 4 shows the broad-band X-ray luminosity 
against the ^T-band absolute magnitude of the X-ray sources 
associated to cluster members, colored according to whether 
they are detected at 24/im (magenta) or not (black) points. X- 
ray point sources associated to field galaxies at 0.15<z<0.30 
are shown in blue. The predicted relations between L x and 
M^ for LMXBs (dot-dot-dot-dashed line; Kim & Fabbiano 
2004) and thermal emission from the diffuse gas halo (dot- 
dashed line; Sun et al. 2007) are both shown, and both lie con- 
siderably below any of the X-ray point sources. From these 
analyses it seems likely that all of our X-ray point sources 
are primarily powered by central AGN rather than any other 
source. 

Examining Fig. 4 it seems that the most X-ray luminous 
AGN are infrared-bright, while many of the lower luminos- 
ity X-ray AGN are massive galaxies not detected at 24/im. 



A Wilcoxon-Mann- Whitney U-test finds that those 24 X-ray 
AGN with SFRs>2M yr _1 have systematically higher L x 
than those with lower SFRs (including 24/im non-detections) 
at the 2.2-cr level (/V=0.014). 

We find no evidence for evolution within the LoCuSS sam- 
ple. Splitting the cluster sample into two according to red- 
shift, and considering only X-ray AGN with L x >10 42 erg s -1 
where we should be reasonably complete (barring the clus- 
ter cores), we find 14 X-ray AGN in the 13 z<0.227 clusters 
(f AGN =0.1%^ x %) and 14 X-ray AGN in the 13 z>0.227 clus- 
ters (fAGN=0-S9^l%). Here the AGN fractions are derived 
for the M^<-22.6 cluster galaxies covered by our Chandra 
data. We do find marginal evidence for a correlation with the 
(thermo)dynamical state of the cluster, as described by a, the 
slope of the logarithmic gas density profile at 0.04 rsoo, or the 
cool-core strength (Sanderson et al. 2009; Smith et al. 2010a). 
Splitting the cluster sample into two according to a, we find 
20 X-ray AGN in the 14 cool core clusters with a< - 0.5 
(f AGN =l.l3^ 3 2 l%), and just 8 X-ray AGN in the 12 non-cool 
core clusters (a>-0.5 (f AGN =0.51^l%). Thus X-ray AGN 
are more frequent in clusters with strong cool cores at the 2 a 
level than those with weak or no cool cores. This higher fre- 
quency of X-ray AGN in cool core clusters is predominately 
due to those with little or no star formation (SFRs<2M yr _1 ): 
9 such X-ray AGN are found in the 14 cool core clusters 
(/agn, low 5F/e=0.51^7%) but just 2 in the 12 non-cool core 
clusters {/agn, low sfr=0.11^%). 

4.1. The X-ray luminosity function of cluster AGN 

Figure 5 shows the broadband (0.3-7 keV) X-ray luminos- 
ity function (LF) of X-ray AGN belonging to the 26 clusters. 
This confirms that ~ 1-2 AGN with L x > 10 42 erg s -1 are found 
per cluster, consistent with Gilmour et al. (2009). Of these, 
27/28 are moderate-luminosity AGN below L^^IO 44 ergs -1 
(Hasinger et al. 2005). We find just one X-ray AGN with 
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Fig. 6. — Stacked observed phase-space diagram (r pro j/rsoo vs (vi os -<v>)/cr v ) of galaxies combining all 26 clusters in our sample. Each gray solid point 
indicates a spectroscopic cluster member, while open points indicate field galaxies. The red curves indicate the \-o velocity dispersion profile in radial bins 
containing 200 galaxies. X-ray AGN are indicated by larger magenta (black) symbols according to whether they are detected (not-detected) with Spitzer at 24/im. 
The size of the magenta/black symbols correspond to X-ray luminosities, Lx<l0 42 ergs -1 (smaller), 10 42 <Lx<10 43 ergs -1 (mid-sized), and Lx>10 43 ergs -1 
(largest). Blue symbols indicate sources spectroscopically classified as broad-line QSOs. The open magenta symbols and error bars indicate the velocity 
dispersion of the cluster X-ray AGN. The dot-dot-dashed box indicates the region with r<0.4r5oo and |v/ O5 -<v>|/a v <0.80 devoid of X-ray AGN. 



quasar-like luminosities (L X >10 erg s ) in our sample, im- 
plying that they are extremely rare within low redshift clusters 
(<0. 13 per cluster; Gehrels 1986). Note that none of the 32 X- 
ray sources lacking redshifts would have Lx>2x 10 43 erg s -1 
if placed at the respective cluster redshifts - i.e. the brightest 
two luminosity bins are robust to spectroscopic incomplete- 
ness. As discussed in § 2.2 not all of the clusters have Chan- 
dra imaging sensitive to X-ray AGN with Lx=10 4L5 ergs -1 . 
We account for this by weighting each X-ray AGN in the 
lowest-luminosity bin (41.5<logLx(ergs -1 )<42.0) by the 
number of clusters in which it could be detected. 

The luminosity function of X-ray AGN in clusters has the 
same shape as, and a normalization a factor ~7.5x higher 
than, that of field X-ray AGN (dashed curve in Fig. 5). The 
field luminosity function is adapted from Hasinger et al. 
(2005). First we multiplied their luminosity function by the 
volume probed by our survey (3.1 x 10 4 Mpc 3 ). Note that for 
a typical cluster in our survey, the depth associated with the 



redshift range over which cluster membership, and thus clus- 
ter volume, is defined (e.g. z=0.195-0.222 for A 209) is of 
order &~0.02 or ~90 Mpc, i.e. much larger than the physical 
sizes of the clusters themselves (r2oo~ 1.5-2.5 Mpc). Second 
we took the analytical fit from their soft (0.5-2.0 keV) X-ray 
LF for 0.20<z<0.40 AGN, evolved it to our mean cluster red- 
shift (z=0.225) based on their best-fit evolution in the space 
density of AGN, (1+z) 4 ' 90 , then corrected the X-ray luminosi- 
ties from their soft band to our broad band assuming a T=1.7 
power-law spectrum. 

The higher normalization of the cluster AGN luminosity 
function is unsurprising because clusters are overdense re- 
gions of the universe. To account for this, and thus achieve 
a more direct comparison between cluster and field AGN 
populations, we use the data from our spectroscopic survey 
to estimate that the redshift- space density of cluster galaxies 
brighter than M K =-23.1 is a factor of ~23 x higher than that 
of field galaxies, when measured over the regions observed 
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with Chandra. The net normalization of the cluster AGN lu- 
minosity function after accounting for the overdense nature 
of the clusters is therefore a factor of ~3 lower than the field 
AGN luminosity function. 

4.2. Dynamical analysis of the X-ray AGN 

In Figure 6 we show the observed stacked caustic diagram 
of all 26 clusters, in which the cluster-centric radii are normal- 
ized by the r^oo determined from the Chandra data, and the 
LOS velocities are scaled by the velocity dispersion (cr v ) of all 
cluster members within r 2 oo (taken to be 1 .5 r^; Sanderson & 
Ponman 2003). The larger symbols indicate the X-ray AGN 
colored magenta (black) according to whether they are (not) 
detected in the mid-infrared with Spitzer, while the size of the 
symbols indicate the X-ray luminosity. The small gray points 
indicate the remaining cluster galaxies, forming the character- 
istic trumpet- shaped caustic profile. 

It is immediately apparent that the X-ray AGN do not trace 
the same distribution in the caustic diagram as the remain- 
der of the cluster galaxy population, but preferentially trace 
the caustics, suggestive of an infalling population (§ 3). The 
24/im-detected AGN trace the caustics better than the non- 
24/im detections (we will return to this later). 

The X-ray AGN also notably avoid the area with low 
cluster-centric radii and relative LOS velocities, although 
we note that the X-ray emission from the ICM affects our 
ability to detect X-ray AGN at these low radii. Other 
spectroscopically-classified broad line QSOs (blue symbols) 
also appear to avoid these central regions. There are no X-ray 
AGN with r<0.4r 5 oo and |v/ os -<v>|/cr v <0.80 (within the 
dot-dot-dashed box in Fig. 6). In this same region of phase- 
space we have 423 cluster galaxies with M^<-23. 1, of which 
we may expect 6. 1 ±2.5 to be X-ray AGN by chance, based on 
the overall detection rate of X-ray sources (i.e. not just those 
with L X >10 42 ergs -1 ) among the general cluster population 
(with Chandra coverage). The probability of detecting none 
is just 0.002. This region of phase-space should be dominated 
by those galaxies accreted earliest into the clusters (§3). 

We do note a couple of biases which may contribute to 
this apparent void. Firstly, in the outer regions (r>0.4r5oo) 
the luminosities of X-ray AGN with |v-<v>|/a v <0.80 are 
~3 x lower than their counterparts with higher velocity off- 
sets. Secondly, the ability to detect faint X-ray AGN using 
wavdetect is reduced in the cluster cores due to increased 
photon noise from the ICM (Fig. 1). Hence if this trend for 
fainter X-ray sources to lie predominately at lower velocity 
offsets seen at r>0.4r5oo were to hold also in the core region, 
we may well be systematically missing them. However, there 
should be no observational bias with respect to LOS velocity. 
Indeed there are many X-ray AGN detected near the cluster 
cores (r<0.4r5oo), but all have large LOS velocities (~1 000- 
3 700 km s -1 ) with respect to that of the cluster. Moreover, 
the fact that the optically- selected AGN also avoid the cluster 
core (blue symbols in Fig. 6), yet shouldn't be affected by the 
same radial bias as the X-ray AGN, supports the view that this 
void is real. 

The mean velocity dispersion of the X-ray AGN is 
1.51 ±0.11 a v (shown by open magenta symbol with error 
bars), significantly higher than that observed for the gen- 
eral cluster population (red curve). Using the Wilcoxon- 
Mann- Whitney non-parametric U test, we find that the ab- 
solute LOS normalized velocities \vi os -<v>\/a v of the 48 
X-ray AGN are systematically higher than the remainder of 
the cluster population (just those covered by our Chandra 




line-of-sight velocity / a v 



Fig. 7. — The stacked velocity distributions (normalized by a v ) of X- 
ray AGN (thick magenta curve) and inactive cluster members (red solid his- 
togram) within r2oo- The blue dot-dashed histogram shows the stacked ve- 
locity distribution of simulated galaxies on their first infall into the 30 most 
massive clusters from the Millennium simulation that lie within a projected 
radius of r2oo and a physical (3D) cluster-centric distance of 2 r2oo- 



images) at the 4.66-a significance level (P[/~1.6x 10~ 6 ). In 
contrast we find no statistical difference in the radial distri- 
bution of the X-ray AGN from the remainder of the clus- 
ter population. Comparing visually the X-ray AGN detected 
with Spitzer (magenta points) and those not detected (black 
points), there is a suggestion that the Spitzer detections have 
systematically higher absolute line-of-sight velocities. A 
Wilcoxon-Mann- Whitney test finds those 23 X-ray AGN with 
obscured SFRs>2M© yr _1 (or L TIR >2x 1O 1O L ) have sys- 
tematically higher |v-<v>|/a v than those with lower SFRs 
at the 2.3a level. Similarly, the more X-ray luminous AGN 
(indicated by the larger symbols) appear to have higher ve- 
locity offsets than their lower luminosity counterparts. Again, 
a Wilcoxon-Mann- Whitney test finds that the 29 X-ray AGN 
with L x > 10 42 erg s -1 have systematically higher | v-<v> | / cr v 
than those with lower X-ray luminosities at the 2.5a level 
(/V=0.005). 

The velocity distributions of X-ray AGN and inactive (not 
X-ray point sources or 24/im-detected) cluster galaxies for 
our stacked cluster sample are compared directly in Figure 7 . 
While the velocity distribution of inactive galaxies within r2oo 
(solid red histogram) can be described approximately as a 
Gaussian function of width a v , the velocity distribution of X- 
ray AGN (magenta curve) appears more consistent with a flat, 
top-hat profile than a Gaussian. The kurtosis of the velocity 
distribution of X-ray AGN (1.965) is significantly lower than 
that expected for a Gaussian distribution (3.0) at the 98.9% 
confidence level, but close to expectations from a uniform dis- 
tribution (1.8). 

We may be concerned that many of the X-ray AGN at high 
absolute LOS velocities are not in fact cluster galaxies, but 
interlopers from the general field population. However, we 
identify just 14 X-ray AGN in the general field population at 
0.15<z<0.30 in the same Chandra images, that is after ex- 
cluding the redshift ranges deemed to contain the cluster pop- 
ulations. Based on this field X-ray AGN space density, we 
would expect ~3 field X-ray AGN to be interlopers within the 
cluster population by chance for the whole survey of 26 clus- 
ters, giving an estimated field contamination of 6%. Hence it 
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Fig. 8. — Color-magnitude (g-i/Mjc) diagram of cluster galaxies (small symbols), color coded according to their specific-SFRs (combining the 24 ^mi-based 
SFRs estimated using the calibration of Rieke et al. (2009) with stellar masses derived from M^). The g—i color of each cluster galaxy is k-corrected to z=0.2 
enabling all of the clusters to be stacked onto a single combination of a red sequence of passive galaxies not detected at 24/mi (red points) and blue cloud of 
star-forming galaxies (green/blue points). X-ray AGN are indicated by larger magenta (black) symbols according to whether they are detected (not-detected) 
with Spitzer at 24/zm. Large blue symbols indicate sources spectroscopically identified as broad-line QSOs. Seven X-ray AGN do not have SDSS coverage and 
three more are fainter than M £+2 and are not plotted. 



seems unlikely that the bulk of the X-ray AGN located along 
the caustics in Fig. 6 are interlopers and there by chance due 
to projection effects, but are instead a genuine population of 
galaxies inf ailing into the clusters for the first time. The high- 
est velocities of cluster galaxies occur for those galaxies on 
their first cluster infall, having passed through the virial ra- 
dius on predominately radial orbits. 

From this comparison to the field population, and the ob- 
served X-ray luminosity function (§ 4.1), the density of X- 
ray AGN in redshift- space within clusters appears ~8-16x 
higher than that for the general field. This does not mean that 
the encountering of the cluster environment is somehow "trig- 
gering" the X-ray AGN, but simply reflects the much higher 
space densities of galaxies as a whole within clusters. Indeed, 
the fraction of cluster members with M#<-23.1 identified as 
an X-ray point source (of any L x ) is 1.4±0.2%, is marginally 
lower than that observed for the 0.15<z<0.30 field galaxy 
population in the same Chandra images (1.8±0.5%). 

4.3. Optical colors, spectra and morphologies 

We have spectral coverage for 42 of the 48 X-ray AGN from 
our ACReS MMT/Hectospec dataset, the remaining six red- 
shifts coming from the literature. Of those for which we have 
spectra: 16 (38.1^g \%) show broad-line Ha emission of Type 
I Seyferts/QSOs; 17 more (40.5^%) show narrow-line Ha 



emission which may be due to either nuclear activity or star 
formation; while 9 (21.4^j%) appear passive, lacking any 
clear emission lines. All of the Type I AGN are detected at 
24/im as opposed to just 3/9 of the passive galaxies. 

Figure 8 shows the optical color-magnitude (g-i/M K ) dia- 
gram of cluster galaxies for the 22 clusters for which SDSS 
DR7 ugriz photometry is available. The IR-dim X-ray AGN 
(black symbols) lie on or close to the cluster red sequence, 
suggesting little impact from the nuclear activity on galaxy 
color. The 24/im-detected X-ray AGN (magenta symbols) 
mostly lie within the blue cloud, suggesting star-formation 
concurrent with nuclear activity. These objects appear blue 
due to their continuum flux, rather than AGN emission lines 
affecting their colors. 

The morphological classification correlates strongly with 
location within or below the cluster red sequence. All but 
one of the 24/im-detected X-ray AGN below the red sequence 
are spirals or unresolved, while all of the >L* X-ray AGN 
along red sequence are E/SOs. Overall, we estimate that 24/48 
(50±8.2%) of the X-ray AGN are early-types (E/SOs), 8 are 
undisturbed spirals (16.7^%), 9 disturbed spirals/mergers 
(18.8!^%) and 7 unresolved objects (14.5!™%), which 
are either low-luminosity galaxies or QSOs. To determine 
morphological classifications we have high-resolution optical 
imaging from the Hubble Space Telescope for 14 out of the 
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Fig. 9. — Thumbnail HST images of all 14 X-ray AGN from our cluster sample for which HST imaging is available. The visual morphology and spectral 
classification are indicated for each host galaxy where available. 



48 X-ray AGN identified as cluster members, thumbnails of 
which are shown in Fig. 9. Eleven of the 14 are detected with 
Spitzer, five of which appear to be spirals, plus one merger. 
The three not detected with Spitzer appear as early-types (two 
SOs, one E). For the remaining 34 X-ray AGN we have high- 
quality Subaru optical imaging which affords visual morpho- 
logical classifications, albeit with difficultly in distinguishing 
SOs and early-type spirals as spiral arms get smoothed out, so 
some spirals may be mis-classified as early-types. The dis- 
tribution of morphologies shown in the HST thumbnails are 
thus representative of the whole sample. 

All of the undisturbed and disturbed spirals and mergers are 
detected at 24/im, as are 5/7 of the unresolved objects. Of the 
32 X-ray AGN detected at 24/im, 17 are spirals (both undis- 
turbed/disturbed), five are unresolved and the remaining 10 
are early-type galaxies, of which all five with HST imaging 
were classified as SOs. We note that four of these 10 early- 
types show broad-line optical emission. Of those 16 X-ray 
AGN not detected with Spitzer, 14 are morphologically iden- 
tified with early-types, the other two being unresolved. 

5. DISCUSSION 

We have studied the distribution of X-ray AGN in 26 
massive clusters at 0.15<z<0.30 from LoCuSS, combining 
Chandra imaging sensitive to X-ray sources of luminosity 
L X ^10 42 ergs -1 at the cluster redshift, with extensive and 
highly complete spectroscopy of cluster members. In total 
48 X-ray AGN were identified among the cluster members, 
with luminosities 2 x 10 41 - 1 x 10 44 erg s -1 . 

5.1. The dynamical status of X-ray AGN within clusters 

The principal result of our analysis of the X-ray AGN pop- 
ulation within the 26 clusters was shown in Fig. 6, which 
showed their location within the stacked caustic diagram, 
(vios~<v>) / o- v versus r pro j/rsoo- They are preferentially lo- 
cated along the caustics, suggestive of an infalling population. 
They also appear to avoid the region with the lowest radii and 
relative velocities (r<0.4r5oo; |v/ OJ -<v>|/cr v <0.80; dot-dot- 
dashed box in Fig. 6), which is dominated by the virialized 
population of galaxies accreted earliest by the clusters. This 
dynamical behavior also shown by the optically- selected Type 
I Seyferts/QSOs, which also make up 38.1+|;|% of the X-ray 
AGN sample. Based on the simulated stacked caustic plots 
of the 30 most massive clusters in the Millennium simula- 
tions at z=0.21 (Fig. 3), 86 (71) percent of all galaxies with 
r<0.4r5oo and |v/ os -<v>|/cr v <0.80 were accreted more than 
1.5 (3.7) Gyr previously, while just 5% are infalling galaxies 
yet to pass within rsoo- 



The velocity dispersion of the X-ray AGN population is 
~50% higher than the overall cluster population. This kine- 
matic segregation is significant at the 4.66-a level. Consid- 
ering a simple kinematical treatment of infalling and virial- 
ized cluster galaxies in a cluster- scale potential well leads to 
| T/V | « 1 for infalling galaxies and | T /V | ~ 1 /2 for the virial- 
ized population, where T and V are the kinetic and potential 
energies. Thus, the velocity dispersions of infalling and virial- 
ized galaxies are related by (Ji n faii~^c>viriai (Colless & Dunn 
1996), close to the ratio observed between the X-ray AGN and 
inactive populations. For our simulated cluster galaxies, the 
highest velocity dispersions are seen for galaxies which have 
been accreted by the cluster within the last 700 Myr, passing 
through r5oo for the first time between 0.21<z<0.28 (light 
blue/turquoise points in Fig. 3), reaching ~1.5a v , i.e. sim- 
ilar to that seen for our X-ray AGN. The velocity dispersion 
rapidly drops back to < 1 .0a v for galaxies accreted earlier into 
the cluster, and is also too low (~1.0cr v ) for those yet to pass 
within r5oo for the first time. 

The velocity histogram (Fig. 7) shows a relatively flat distri- 
bution, unlike the approximately Gaussian profile of the inac- 
tive galaxy population, as confirmed by its measured kurtosis 
being lower than that expected for a Gaussian distribution at 
the 98.9% confidence level. Sanchis, Lokas & Mamon (2004) 
show that high velocity dispersions and low kurtosis values 
are predictors of galaxies on predominately radial orbits, at 
least in the cluster cores. This flat- topped distribution is also 
seen for galaxies within a projected radius of r2oo that are in- 
falling into our simulated clusters for the first time (blue dot- 
dashed histogram in Fig. 7), and yet to reach the pericenter of 
their orbit. 

Both these velocity profiles and the kinematic segregation 
of the X-ray AGN in the stacked caustic profiles are similar to 
those observed for spiral and/or emission-line galaxies within 
clusters (e.g. Biviano et al. 1997, 2002; Boselli & Gavazzi 
2006; Haines et al. 2010). From the first dynamical studies 
of cluster galaxies, the velocity dispersions of spiral galaxies 
have been found to be systematically higher than early-types 
(Tammann 1972; Moss & Dickens 1977). Based on much 
larger cluster samples, Biviano et al. (1997), Adami, Biviano 
& Mazure (1998) and Aguerri et al. (2007) all found that 
the stacked velocity dispersions of blue/emission-line galax- 
ies to be on average 20% higher than the remaining inactive 
galaxies. Biviano et al. (1997) also examined the kinemat- 
ics of spectroscopic AGN in local clusters, finding them to 
have systematically higher velocity offsets than the inactive 
galaxy population. Biviano et al. (2004) later showed that if 
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Fig. 10. — Comparison of the locations in the phase-space diagram {r pro j/r^Q vs (yi os -<v>)/ cr v ) of infalling (left panel) and back-splash (right panel) cluster 
galaxies from the stacked cluster sample from the Millennium simulation. Infalling galaxies are defined as those with negative radial velocities which have yet 
to reach pericenter in their orbit around the cluster. They are color coded according to their physical (3D) distance from the cluster center from red (r<0.25 r^oo) 
to blue (r>2.0r5oo). Back-splash galaxies are defined as those which have passed through the pericenter for the first time, and their orbits are now taking them 
back out from the cluster towards apocenter. Their symbols are instead color coded as in Fig. 3. 



early-type galaxies are assumed to have isotropic orbits within 
clusters as supported by the Gaussian shape of their velocity 
distribution, the kinematic properties of late-type spirals are 
inconsistent with being isotropic at the >99% level. Instead 
they indicate that spirals and emission-line galaxies follow ra- 
dial orbits in clusters, pointing towards many of them being 
on their first cluster infall. 

The cluster environment is known to strongly impact the 
evolution of member galaxies, transforming infalling star- 
forming spiral galaxies into passive early-type galaxies, as 
manifest by the SF-density (Dressier, Thompson & Schect- 
man 1985) and morphology-density (Dressier 1980) rela- 
tions. The kinematic segregation of star-forming and pas- 
sive galaxies in clusters is a key empirical foundation for our 
understanding of this process, by revealing that, while the 
passive galaxies are a virialized population that have resided 
within the cluster for many Gyr, the late-type star-forming 
galaxies are much more recent arrivals, and indeed many have 
yet to encounter the dense intra-cluster medium. It is be- 
lieved that this arrival of the spiral galaxy within the cluster 
and its passage through the ICM removes its gas supply, lead- 
ing to the subsequent quenching of its star formation (Boselli 
& Gavazzi 2006, and references therein). That the kinematic 
signature of the X-ray AGN within clusters is so similar to 
that of star-forming spirals, indicates that they are also recent 
arrivals into the cluster. Their apparent absence in the cen- 
tral regions dominated by those galaxies accreted earliest into 
clusters, confirms that the harsh cluster environment strongly 
suppresses radiatively-efficient nuclear activity. 

We find further dynamical support for the ongoing sup- 
pression of nuclear activity in galaxies as they arrive in the 
clusters, with both X-ray bright (L^>10 42 erg s -1 ) and IR- 
bright (L77tf>2xlO lo L ) sub-samples of X-ray AGN show- 
ing higher velocity dispersions than their X-ray dim and IR- 
dim counterparts at the >2a level. This is consistent with the 
nuclear activity responsible for the X-ray and infrared emis- 
sion being gradually shut down as the host galaxies are ac- 
creted into the cluster, losing their orbital kinetic energy as 
part of the long dynamical process of becoming virialized 
cluster members. We note however, that the typical duty cycle 
of AGN of 10 7 -10 8 yr is much lower than the >l-2x 10 9 yr 
cluster crossing time-scales which are likely required for such 



velocity segregation to occur (e.g. Gill, Knebe & Gibson 
2005). Instead we suggest that the velocity segregation re- 
flects the slow decline in the gas contents of the host galaxies 
available for intermittently fuelling nuclear activity. 

We attempt to confirm this viewpoint by going back to the 
stacked caustic plot for galaxies orbiting the 30 most massive 
clusters in the Millennium simulation discussed in § 3. In 
Figure 10 we replot the caustic diagram of Figure 3, but now 
split the galaxies into specific sub-populations based on their 
orbital paths. 

In the left-hand panel we show only those galaxies which 
are on their first infall into the cluster, that is they have neg- 
ative radial velocities and are yet to reach pericenter for the 
first time. Each point is color coded according to its physical 
(3D) distance from the cluster center in units of rsoo- We see 
that these galaxies show a wide range of line-of- sight veloci- 
ties, evenly distributed over the phase-space with a slight pref- 
erence towards the caustics at ~±2a v , as shown previously 
by their flat- topped velocity dispersion in Figure 7. This was 
found to be the sub-population with the distribution of galax- 
ies in the caustic diagram which was the most qualitatively 
similar to that seen for the X-ray AGN detected by Spitzer 
(the magenta points in Fig. 6). We note that although there is 
no concentration of galaxies towards low relative LOS veloc- 
ities and radii (as seen for the earliest accreted galaxies), they 
do not entirely avoid this region either (as our observed X-ray 
AGN apparently do). This could only be achieved by remov- 
ing those galaxies within physical radii of <0.5 rsoo (red, or- 
ange points) and also those at >2r5oo (blue points). However 
this would also impact the required presence of galaxies at 
±2a v and low projected cluster-centric radii (<0.4r5oo). We 
believe that many of the X-ray AGN are near pericenter as 
required to achieve their high LOS velocities. 

The galaxies with the highest LOS velocities along the 
caustics at ~±2a v are found to be infalling galaxies at ~0.25- 
1-0 r5oo (orange, yellow points). Such high LOS velocities are 
notably not seen for those infalling galaxies at large physical 
distances from the cluster (>2r5oo; blue points). These galax- 
ies have detached from the Hubble flow but haven't fallen suf- 
ficiently far into the cluster gravitational potential well to at- 
tain the very high infall velocities required to be located along 
the caustics, resulting in them having relative line-of-sight ve- 
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locities within ~1.0cr v . This could suggest that the nuclear 
activity is being triggered in these infra-red bright galaxies as 
they are accreted into the cluster for the first time, at ~ 1-2 rsoo 
from the cluster center, when the galaxies have reached close 
to their maximal velocities, or indeed at pericenter passage. 
This phase could be induced by galaxy harassment whereby 
inf ailing spiral galaxies undergo frequent high-speed fly-by 
encounters with other cluster members, or tidal shocks as they 
pass through the cluster core, driving instabilities that funnel 
gas into the central regions triggering nuclear activity (Moore 
et al. 1996). This would be consistent with our finding that all 
but one of the IR-bright X-ray AGN are spirals or unresolved. 
Unfortunately, as our Chandra data do not extend into the in- 
fall regions beyond 1-1.5 rsoo, we cannot easily test whether 
there are extra X-ray AGN among the infalling population just 
at these cluster-centric radii (as opposed to further out in the 
inf all regions). 

In the right-hand panel we show the back- splash popula- 
tion of galaxies which have passed through the cluster core 
for the first time and are on their way back out of the clus- 
ter. That is they are between first pericenter and apocenter 
along their orbits and have positive radial velocities. In terms 
of cluster-centric radius this population peaks at ~1 rsoo, be- 
coming steadily less frequent towards the lowest projected 
radii. This is due to them having their highest orbital veloc- 
ities at pericenter, slowing down as they approach apocen- 
ter and hence spending more time near this point. Although 
the spread of LOS velocities for the back-splash population is 
similar to that for the infalling galaxies, they are more concen- 
trated towards low LOS velocities, as seen previously by Gill, 
Knebe & Gibson (2005) and Mahajan et al. (2011). We be- 
lieve this distribution more closely resembles that seen for our 
X-ray AGN not detected with Spitzer (black points in Fig. 6). 
This back-splash population is observed ~ 1 Gyr later after ac- 
cretion into the cluster than the infalling population, as in- 
dicated by their predominately green symbols in the Figure. 
This suggests that as the galaxy passes through the cluster 
core for the first time, the gas available for star formation and 
its associated infrared emission is removed via ram-pressure 
stripping, while the gas available for nuclear activity, being 
deeper in the galaxy potential is harder to sweep out, allow- 
ing nuclear activity to continue past pericenter passage. The 
different morphologies of the galaxies hosting the IR-bright 
(mostly spirals) and IR-dim (mostly E/SOs) X-ray AGN sug- 
gests however that the primary factor in determining whether 
the host galaxy is IR-luminous is the global availability of gas 
and dust in the galaxy host itself, IR-dim X-ray AGN being 
hosted in early-type galaxies which have sufficient gas in the 
nucleus to fuel the AGN, but not enough on larger scales to 
feed obscured star formation. 

The results of our dynamical analysis of X-ray AGN are in 
marked contrast to the comparable survey of (Martini et al. 
2007) who found that the spatial and kinematical distribution 
of their X-ray AGN sample were fully consistent with being 
drawn from the inactive cluster population. Given that their 
sample size is similar to ours, it is surprising that the results 
are so different and apparently inconsistent. The primary dif- 
ference between the surveys is that their X-ray AGN are sys- 
tematically less luminous than ours. Hence, in the scenario 
described above in which nuclear activity is slowly quenched 
when the host galaxy is accreted into the cluster, the low- 
luminosity X-ray AGN that dominate the sample of Martini 
et al. (2007) may have witnessed the cluster environment for 
longer than those in our sample, and so are less kinematically 



distinct from the general cluster population. Interestingly, 
Ajello et al. (2012) find that the fraction of Seyfert 2 objects 
among their all- sky hard X-ray AGN sample is much higher 
within the largest concentrations of matter in the local Uni- 
verse, suggesting that the broad-line Type I AGN are being 
preferentially shut down in these dense environments, consis- 
tent broadly with our kinematic segregation of IR-bright/IR- 
dim X-ray AGN. 

We note finally that despite our result that the X-ray AGN 
are an infalling population (or just after pericenter) which is 
strongly suppressed by the cluster environment, they repre- 
sent a ~8-16x overdensity in redshift space with respect to 
the general field, as revealed by the X-ray luminosity func- 
tion of cluster AGN (§4.1). This is fully consistent with the 
statistical overdensities of X-ray AGN previously observed 
toward clusters with respect to non-cluster fields (e.g. Cappi 
et al. 2001; Molnar et al. 2002; Ruderman & Ebeling 2005). 
This does not imply that they are somehow triggered by the 
cluster environment, but instead reflects the overall significant 
(~ 1 0-200 x) overdensities of galaxies in the inf all regions of 
clusters, that have yet to encounter the ICM. Indeed the red- 
shift space density of normal cluster galaxies in our sample 
is ~23x higher than of field galaxies, higher than the re- 
spective overdensity of cluster X-ray AGN. It is still possible 
that there is an increased frequency of nuclear activity among 
these infalling galaxies, due to encounters between galaxies 
in the connecting filamentary web, or pre-processing within 
infalling groups, although our lack of Chandra coverage of 
these inf all regions prevents us from measuring this. Nuclear 
activity may also be triggered by their interaction with the 
cluster itself, either as they pass through virial shocks, via 
compression of gas onto the nucleus in the early stages of 
ram-pressure stripping, or tidal shocks as they pass through 
cluster pericenter. 

5.2. Comparison to results from clustering analyses 

A complementary approach to constraining the typical en- 
vironment of AGN, and gain insights into the physical condi- 
tions of the accretion onto SMBHs, is to measure the spatial 
clustering of X-ray AGN to estimate the average DM halo 
mass harboring the AGN and the relative frequency of being 
hosted by central or satellite galaxies. Numerous clustering 
analyses have revealed a consistent view that X-ray AGN are 
hosted in DM halos of masses ~2x 1O 13 M , typical of poor 
groups, at all redshifts up to z~2 (Hickox et al. 2009; Cap- 
pelluti et al. 2010; Allevato et al. 2011; Miyaji et al. 2011). 
Starikova et al. (2011) show that they are predominately lo- 
cated at the centers of DM halos of mass >6x 1O 12 M , and 
tend to avoid satellite galaxies in comparable or more mas- 
sive halos, fixing the limit to the fraction of AGN in satel- 
lite galaxies to be <12% (90% confidence limit). Similarly, 
Hickox et al. (2009) find the clustering of X-ray AGN to 
be consistent with that of typical galaxies on scales of 1- 
lOMpc, but significantly anti-biased on small scales (0.3- 
1 Mpc), which they explain as due to their preferential loca- 
tion within central galaxies. Tasse, Rottering & Best (2011) 
also find that on 450 kpc scales the most luminous X-ray AGN 
(L^>10 43 ergs -1 ) are found in underdense environments in 
comparison to normal galaxies of the same stellar mass. 

In direct contrast, Mountrichas & Georgakakis (2012) 
found X-ray AGN to be more clustered than galaxies at 
all scales, with no evidence for anti-bias on small scales. 
While they also obtained a typical host DM halo mass of 
~lO 13 /z _1 M , the stellar masses of their host galaxies were 
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lower than expected for the typical central galaxy of such a 
halo, suggesting that they are associated to satellite galaxies. 
Miyaji et al. (2011) found excess clustering on small scales 
between X-ray AGN and luminous red galaxies, consistent 
with a significant fraction of X-ray AGN being hosted in satel- 
lite galaxies. 

Miyaji et al. (2011) also indicate that the AGN frac- 
tion in groups and clusters declines with halo mass. 
Arnold et al. (2009) also find that the X-ray AGN frac- 
tion (L^>10 41 ergs -1 ) declines by a factor two from groups 
(fAGN=0mi+™f 4 ) to clusters (/ag^O.047^) from their 
study of 16 systems at 0.02<z<0.06 . This trend remained 
when considering only early-type galaxies, and so was inde- 
pendent of any differences in the morphological mix between 
groups and clusters. Hwang et al. (2012) also observed a de- 
cline in AGN fraction from groups to clusters (including at 
fixed morphology) for optically- selected AGN. 

The observed dynamics of the X-ray AGN in our cluster 
sample are consistent with their being entirely drawn from 
an infalling population, with essentially no component com- 
ing from virialized satellite galaxies within the cluster halo. 
These results are consistent with the clustering analyses of 
Hickox et al. (2009) and Starikova et al. (2011), associating 
the X-ray AGN with central galaxies. Their reported anti-bias 
on small scales then reflects our observed suppression of nu- 
clear activity among the virialized cluster galaxy population. 
It is instead hard to reconcile our results with the findings of 
Miyaji et al. (2011) and Mountrichas & Georgakakis (2012) 
that X-ray AGN being even more clustered on small scales 
than the already strongly clustered luminous red galaxies, as 
in this case they should be preferentially located in the cluster 
cores rather than avoiding them as we find. 

5.3. The properties of the galaxies hosting X-ray AGN in 

clusters 

Studying the physical properties of the galaxies which host 
AGN provide complementary clues as to their accretion pro- 
cesses. The stellar masses, bulge masses or stellar velocity 
dispersions can be used to estimate the mass of the central 
SMBH via the known tight correlations, from which the Ed- 
dington ratio can be derived. In the absence of direct measures 
of the Hi or molecular gas contents, galaxy colors can be used 
as a proxy for star-formation history, and hence the availabil- 
ity of gas. Finally, galaxy morphologies provide constraints 
on the stage of any ongoing merging event (plus mass ratio) 
or the presence of bars or disks required for the related secular 
processes. These probes are however more prone to system- 
atics, as the AGN can easily outshine the hosts, affecting not 
only their luminosities but also colors. 

Hickox et al. (2009) found that while radio AGN are hosted 
mainly by massive, red sequence galaxies, the X-ray and 
infrared- selected AGN are instead both found in — L* galax- 
ies, with the X-ray population being preferentially "green val- 
ley" objects, while IR AGN are slightly bluer. Haggard et al. 
(2010) also found that X-ray AGN are much more likely to be 
located within the blue cloud or green valley than the red se- 
quence. Georgakakis et al. (2009) suggest that the color distri- 
bution for X-ray AGN hasn't evolved between z— 0.8 and the 
present day. The optical colors and luminosities of our cluster 
X-ray AGN are fully consistent with these field samples, be- 
ing hosted typically by — L* galaxies located both within the 
red sequence and the blue cloud (Fig. 8). We find marginally 
more X-ray AGN along the cluster red sequence (49 ±9%) 
than Hickox et al. (2009) (39±4%), which could reflect a 



cluster- specific difference in the host properties, such as re- 
duced levels of star formation, but more likely reflects our use 
of g—i color rather than u-r to define the red sequence, and the 
systematically lower X-ray luminosities in our sample. The 
IR-dim X-ray AGN lie on or close to the cluster red sequence 
in Fig. 8, suggesting little impact from the nuclear activity 
on galaxy color. We find the 24/im-detected X-ray AGN to 
lie mostly within the blue cloud, consistent with Hickox et al. 
(2009), and suggesting star-formation concurrent with nuclear 
activity. 

From our spectroscopic analysis of the cluster X-ray AGN, 
we found that 38% show broad-line Ha emission of Type I 
Seyferts/QSOs, 41% show narrow-line Ha emission which 
may be due to either nuclear activity of star formation, while 
21% appear passive, lacking any clear emission lines. These 
fractions are strongly inconsistent with those obtained by 
Martini et al. (2006) from their spectroscopic survey of X- 
ray AGN in clusters. They found that only 4 out of 40 X-ray 
AGN showed clear optical signatures of nuclear activity in the 
form of significant Oil, Olll or broad Ha emission. The re- 
maining 90^% show only modest star formation or appear to 
be passively evolving galaxies. The primary difference be- 
tween the two samples is that the X-ray AGN from Martini et 
al. are systematically — 2.5x less luminous than ours, having 
a mean value of log(Lx)=41 .78 as opposed to our mean value 
of log(Lx)=42.16. We would thus expect the levels of optical 
emission to be correspondingly lower from the unified AGN 
model. Burlon et al. (2011) find that low-luminosity X-ray 
AGN are more likely to be absorbed (N H >10 22 cm -2 ), which 
may further contribute to the lack of optical emission among 
the X-ray AGN sample of Martini et al. (2006). 

Cisternas et al. (2011) found that <15% of X-ray AGN in 
the HST-COSMOS field (with L x -10 435 ergs -1 ; z-0.3-1.0) 
showed any signs of distortions indicative of recent merg- 
ers, and indeed found no statistical difference in the distortion 
fractions between X-ray AGN and inactive galaxies, indicat- 
ing that major mergers are not the most relevant mechanism 
for the triggering of X-ray AGN at z< 1 . Instead they found 
that over 55% of the X-ray AGN are hosted by disk galax- 
ies, and suggest that the bulk of black hole accretion occurs 
through internal secular fuelling processes and minor merg- 
ers. Similarly, Griffith & Stern (2010) found while radio- 
loud AGN are mostly hosted by early-type galaxies, X-ray 
AGN are mostly either disk-dominated (31-46%) or unre- 
solved point sources (31-61%) with few (9-21%) hosted by 
bulge-dominated systems. In the local Universe, Koss et al. 
(2011) found that X-ray AGN are -5-10 x more likely to be 
hosted in spirals (—40%) or mergers (—20%) than inactive 
(non-AGN) galaxies of the same stellar masses. Ellison et 
al. (2011) did show however that some of the nuclear activ- 
ity is triggered by interactions, showing that the AGN frac- 
tion increases by up to —2.5 x for galaxies in close pairs with 
projected separations <10kpc, and that this enhancement in 
nuclear activity is greatest for equal-mass galaxy pairings. 
The morphological composition of our cluster X-ray AGN 
are broadly consistent with these previous studies, with 24/48 
(50%) being early-types (E/SOs), 8 undisturbed spirals (17%), 
9 disturbed spirals/mergers (19%) and 7 unresolved objects 
(14%), which are either low-luminosity galaxies or QSOs. We 
find a marginally higher fraction hosted by early-types, which 
could reflect the general increased prevalence of early-types 
among cluster galaxies, but could also be an overestimate due 
to our inability to robustly distinguish E/SOs and early-type 
spirals from the Subaru imaging. 
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The frequent hosting of X-ray AGN (at least for z<l) in 
otherwise undisturbed massive spirals appears inconsistent 
with the classical association between black hole accretion 
and bulge growth, and indeed the merger paradigm. As such 
we may have expected X-ray AGN to be preferentially hosted 
by early-type galaxies which have the most massive black 
holes, and hence must have had the most nuclear activity in 
the past. Hopkins & Hernquist (2006) have developed a sce- 
nario for the fuelling of AGN in non-interacting spirals/SOs 
in which the accretion of cold gas onto supermassive black 
holes occurs via stochastic collisions with molecular clouds 
(or the inflow of these clouds via disk/bar instabilities). In 
these systems, the only requirement is the availability of cold 
gas within the disk, producing intermittent bursts of accretion 
at high Eddington ratios (~1-10%) and duty cycles of ~1%. 
The resultant feedback from the AGN is expected to have 
negligible impact on the host galaxy's disk and interstellar 
medium. The hosting of low-redshift X-ray AGN within mas- 
sive spirals can thus be simply understood as the requirement 
of an available supply of cold gas. While elliptical galaxies 
host the most massive supermassive black holes, their fre- 
quent lack of cold gas means that there is no fuel for future 
gas accretion or black hole growth, and hence they are less 
likely to be X-ray luminous. Kauffmann et al. (2007) showed 
that early-type galaxies with strong nuclear activity (based on 
their Olll emission) almost always have blue UV-optical col- 
ors and blue outer regions indicative of a reservoir of cold gas 
in the disk capable of feeding the AGN and star formation in 
the outer disk. 



6. SUMMARY 

The key finding from our study of the distribution and host 
properties of 48 X-ray AGN identified from Chandra imaging 
of 26 massive clusters at 0.15<z<0.30 is that they are clearly 
dynamically identified with an inf ailing population. This is 
manifest by their preferential location along the cluster caus- 
tics, complete avoidance of the caustic phase space with low 
relative velocities and cluster-centric radii, and their high ve- 
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locity dispersion and non-Gaussian velocity distribution. The 
optically- selected Type I Seyferts/QSOs in our cluster sam- 
ple show the same kinematical signatures. These provide the 
strongest observational constraints to date that the X-ray AGN 
and optically-selected Type I Seyferts/QSOs found in massive 
clusters are not a virialized population, and few if any can 
have resided within the dense ICM for a significant length of 
time. The cluster environment must very effectively suppress 
radiatively-efficient nuclear activity in its member galaxies. 
The dynamical properties of our X-ray AGN are very similar 
to those previously seen for late-type spiral galaxies in clus- 
ters, and a significant fraction (35. 4+75%) of the galaxies host- 
ing the AGN are morphologically identified as spirals. This 
suggests that they are mostly triggered by secular processes 
such as bar/disk instabilities rather than mergers, although we 
do also find a number of X-ray AGN associated with ongoing 
mergers. 

These results, in conjunction with the previous clustering 
analyses of Hickox et al. (2009) and Starikova et al. (2011), 
all indicate that two key requirements for radiatively-efficient 
nuclear activity to occur in galaxies is that they are the central 
galaxy within their dark matter halo, and have a ready supply 
of cold gas. This latter requirement is supported by the high 
fraction of X-ray AGN hosted by optically blue galaxies, both 
in clusters (51 ±9%) and the general field at z<l. 
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Table 2. Catalog of spectroscopic X-ray AGN identified as cluster members. 



Cluster 


RA 


Declination 


z 




v/a v 


r/r 5 oo 


M K 


SFR 


Morphology 


Spectral 




(J2000) 


(J2000) 




(erg/s) 


(Vega) 


(M /yr) 


Class 


A 1 1 C 

AI I J 


UU.jj.jy.UI j 


+zo.zU.4o.o3 


U.zUzVJ 


I .z4X IU 


Z.ZjjJ 


U.Ujoz 


OA AQA 

-z4.ooU 


U.UUU 




ALrJN. riavJNII 


A115 


00:56:02.271 


+26:27:28.66 


0.19036 


3.15xl0 41 


-0.3121 


0.5686 


-25.832 


0.000 


E/S0 


AGN: Ha;,NII 


A115 


00:56:08.656 


+26:25:10.87 


0.19245 


3.28xl0 41 


0.1115 


0.5996 


-25.317 


0.986 


E/S0 


Passive 


AI I J 


UU.JO.II.4z/ 


+ZO.Z / . j /.Z4 


\jAyyiy 


I. / / X IU 


1 AHH< 

1.4/ /J 


U.oI4j 


-zj.jU4 


j.Ij / 


Unresolved 


rJroadiine i^oU 


A115 


00:56:19.702 


+26:21:53.24 


0.19017 


9.12X10 41 


-0.3506 


0.6700 


-25.895 


1.022 


E/S0 


Passive 


A90Q 


01 -^1 78Q 


1 ^-^9-9^ 1 A 


O 91 8zLfi 


1 « Y 1 042 
1 .Jo X IU 


1 70S6 
1 . /UJO 


1 084^ 
1 .Uo+O 


9zL ^AQ 
-Z4. D^ty 


^ AAl 


Qt->/Tnt9 

op/ini : 


. 

otar torming 


A209 


01:31:50.568 


-13:30:00.38 


0.20255 


6.82xl0 41 


-1.2249 


1.1388 


-22.555 


4.119 


Unresolved 


Broad-line QSO 


A209 


01:32:16.328 


-13:35:38.15 


0.20324 


1.17X10 42 


-1.0978 


0.9475 


-23.104 


4.407 


Merger 


Noisy AGN: Ho;,OII 


A267 


01:52:59.134 


+00:57:04.80 


0.22773 


5.63 xlO 41 


0.0542 


1.2243 


-23.448 


1.702 


Merger 


Noisy AGN: Ha/NII 


AJoJ 


UZ.4o.Z4.o ly 


-U3.JI.4j. JO 


U. Iojj I 


I .UJ X IU 


-I.U4JO 


u.yjoi 


-Z4. /4 / 


1 1 n 
1.11/ 


xi/oU 


Passive 


A383 


02:48:21.892 


-03:34:25.55 


0.1 86 a 


3.16X10 42 


-0.8835 


0.9608 


-25.188 


1.125 


Merger 




A586 


07:32:17.043 


+31:36:50.66 


0.17378 


8.75 xlO 41 


0.8855 


0.1945 


-24.495 


6.798 


Spiral* 


Broadline QSO 


AjoO 


U / .dZ.ZKj.dLZ 


+3I.4I.ZU.OJ 


U.IOJOJ 


/; QQvl fl42 
O.ooX IU 


-I.4UUI 


U.JlOO 


-zj. /Uj 


I.ZJJ 




Passive 


A611 


08:00:56.472 


+36:05:25.02 


0.27960 


9.67 x 10 41 


-1.6040 


0.3828 


-24.929 


4.406 


Sp/Int? 


AGN 


A ^1 1 
AO 1 1 


Uo.UU.jz.yUI 


+JO.UO.ZO.I I 


U.z / /3J 


A QOv 1 cA 1 
O.ooX IU 


-Z.IJ4/ 


n coin 


-ZJ. / ol 


U.UUU 


ii/oU 


Passive 


A611 


08-01 -09 99Q 


-i-^fv09-4S 7^ 


9761 1 


S 97 v 1 042 

J .Z / A IU 


-2 4225 


O 9^84 


9^ ^Hfi 

-ZJ. J JO 


4 7zH 


en* 

OU 




AO 1 1 


OS-01 0£ 1 9S 

UO.Ul .UO. 1ZJ 


+JO.UJ. lO.OJ 


O 977Q0 

U.Z / /7U 


6 67 v 1 041 

O.O / A IU 


9 OOSO 


U. J JDU 


9zL Q78 

-Z+.7 / 


1 8 90S 


opirai 


wll, will 


A611 


08:01:04.616 


+36:02:21.21 


0.27832 


1.14 XlO 42 


-1.9059 


0.3564 


-23.486 


0.000 


SO* 


AGN? Oil, Ha 


A665 


08:29:33.864 


+65:50:59.12 


0.18795 


5.78xl0 41 


1.1041 


1.1408 


-24.614 


0.547 


Spiral 


Just Ha 


A665 


08:30:41.665 


+65:58:10.68 


0.18144 


7.99X10 41 


-0.2650 


0.9509 


-25.987 


0.000 


E/S0 


Broad Ho; 


A665 


08:30:55.923 


+65:58:18.68 


0.18477 


4.90x 10 41 


0.4353 


0.9445 


-25.547 


3.944 


Merger 


Just Ha/NII 


A963 


10:17:00.712 


+39:04:32.76 


0.21039 


7.62 xlO 43 


1.3744 


0.2863 


-25.077 


46.893 


E/S0 


QSO 


A1689 


13:11:16.981 


-01:16:56.63 


0.1898 b 


3.98xl0 41 


0.7731 


0.5845 


-22.457 


0.000 


E/S0 





A1689 


13:11:35.612 


-01:20:12.38 


0.19955 


1.57 xlO 42 


2.3768 


0.1940 


-23.707 


2.727 


Spiral* 


Broadline QSO 


A1689 


13:11:22.161 


-01:23:45.68 


0.1921 b 


1.85 xlO 42 


1.1514 


0.4668 


-22.230 


0.000 


SO* 




A1689 


13:11:45.447 


-01:23:36.07 


0.18685 


1.23 xlO 42 


0.2878 


0.6303 


-24.845 


0.000 


E/S0 


Passive 


A1758 


13:32:44.206 


+50:31:07.82 


0.28733 


1.32x 10 42 


1.3714 


1.3013 


-22.824 


4.088 


Unresolved* 


StarBurst 


A1763 


13:35:42.044 


+41:02:21.79 


0.21816 


5.07 x 10 42 


-2.5409 


0.9325 


-25.177 


8.860 


Spiral 


Broadline QSO 


A1763 


13:35:18.359 


+40:59:50.10 


0.23972 


2.16x 10 42 


1.3333 


0.0235 


-25.524 


7.728 


dusty SO* 


AGN: OIII,NII 


A1835 


14:00:51.309 


+02:59:05.52 


0.25614 


1.99 xlO 42 


0.6792 


1.0289 


-24.176 


4.308 


Unresolved 


QSO 


A1835 


14:01:27.699 


+02:56:06.15 


0.26472 


1.29 xlO 44 


2.0868 


1.0764 


-24.191 


40.856 


Unresolved 


QSO 


A1835 


14-00-48 988 


+09-47-09 59 


95640 


5 39 x 10 41 


0.7218 


0.9855 


-25.236 


0.000 


E/S0 


v J Lcll 1/ LI 1 o L 


A1835 


14:00:54.139 


+02:52:23.38 


0.2617 b 


3.26X10 42 


L5913 


02980 


-22^810 


o!ooo 


Unresolved 




A1835 


14:01:08.836 


+02:44:57.67 


0.2544 b 


5.17X10 41 


0.3937 


1.1788 


-21.318 


0.000 


Unresolved 




A1835 


14:01:16.154 


+02:45:13.22 


0.24468 


1.61 xlO 42 


-1.2009 


1.2298 


-25.501 


0.000 


E/S0 


Passive 


A1914 


14:25:46.628 


+37:49:24.40 


0.17626 


8.85 xlO 41 


2.1850 


0.2861 


-22.887 


0.000 


E/S0 


AGN: Ha/NII 


A2218 


16:35:47.348 


+66:14:44.60 


0.16855 


4.22xl0 41 


-0.9628 


0.2836 


-25.285 


0.000 


E* 


Passive 


A2218 


16-^5-56 090 

iu. jj .yjyyj 


+66-16-15 1^ 

1 VJVJ . 1 VJ . J. J . 1 _> 


0.17932 


1.94x 10 41 


1.2202 


0.5043 


-22.359 


12.967 


SO* 


Nfnkv AGN 


A2218 


16:36:30.913 


+66:15:05.24 


0.18172 


3.75 xlO 42 


1.7067 


0.6773 


-23.671 


1.036 


Spiral* 


Broadline AGN 


A2219 


16:40:05.839 


+46:43:41.43 


0.22891 


3.91 xlO 42 


0.5688 


0.4599 


-23.024 


2.015 


E/S0 


AGN 


A2219 


16:40:09.337 


+46:44:19.55 


0.22993 


1.76 xlO 42 


0.7495 


0.4367 


-24.572 


48.478 


Sp/Int? 


StarBurst 


A2390 


21:53:31.422 


+17:41:33.68 


0.24587 


1.03 xlO 42 


2.9608 


0.1902 


-24.577 


5.658 


Spiral* 


QSO 


A2390 


21:53:25.330 


+17:43:21.75 


0.2242 c 


1.92xl0 42 


-0.8424 


0.4667 


-23.727 


0.000 


E/S0 




A2390 


21:53:37.972 


+17:43:47.46 


0.24071 


6.44 xlO 41 


2.0552 


0.3044 


-23.289 


14.819 


SO* 


QSO 


A2390 


21:53:45.559 


+17:41:47.88 


0.21437 


1.87 xlO 42 


-2.5677 


0.3042 


-25.026 


52.559 


Sp/Int? 


QSO 


A2390 


21:53:45.738 


+17:41:08.46 


0.23483 


6.22xl0 41 


1.0232 


0.3219 


-23.755 


0.000 


E/S0 


Passive 


A2485 


22:48:30.816 


-16:10:30.94 


0.24139 


1.84 xlO 43 


-1.4940 


1.1466 


-25.419 


39.317 


Sp/Int? 


Broadline AGN 


Z1883 


08:42:51.963 


+29:28:25.02 


0.18442 


1.22 xlO 42 


-2.7000 


0.2315 


-23.848 


99.964 


SO* 


QSO 



Note. — Col. 1: Cluster Name. Cols. 2 & 3: The right ascension and declination of the X-ray AGN (J2000). Col. 4: Redshift. Col. 5: Broadband (0.3-7 keV) X-ray luminosity 
of X-ray AGN assuming T=1.7 power law spectrum. Col. 6: Velocity offset relative to the cluster redshift in units of <r v , the velocity dispersion of the cluster. Col. 7: Projected 
cluster-centric radius in units of rsoo- Col. 8: ^T-band absolute magnitude of host galaxy (Vega). Col. 9: SFR of host galaxy in units of M0vr -1 . The SFR is derived from the 24 /im 
luminosity, and is set to zero for Spitzer non-detections. Col. 10: Morphological class of host galaxy. E/S0: Early-type galaxy. Sp/Int?: Spirals with signs of possible interations. 
Starred morphological classifications indicate those host galaxies for which HST imaging was available. Col. 11: Spectroscopic class derived from our MMT/Hectospec data where 
available. 

a Redshift from Rizza et al. (2003). 
b Redshifts from Czoske (2004). 
c Redshift from Crawford et al. (2002). 



